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Chapter  1 
Introduction 


The  study  of  low-frequency  scattering  was  initiated  by  Rayleigh  [1]  in  1897  and  today  the 
term  “Rayleigh  scattering”  is  often  used  instead  of  “low-frequency  scattering”.  Kleinman  [2] 
gives  the  following  definition  of  Rayleigh  scattering: 

IVe  are  dealing  with  Rayleigh  scattering  when  the  far-zone  field  may  be  expanded 
in  a  convergent  series  in  positive  integral  powers  of  the  propagation  constant  k. 

The  starting  point  in  a  Rayleigh  scattering  calculation  is  the  expansion  of  the  unknown 
fields  in  powers  of  k  and  then  the  determination,  from  Maxwell’s  equations  and  boundary 
conditions,  of  the  unknown  expansion  coefficients,  which  are  functions  of  the  geometry  of  the 
scatterer  and  the  angles  of  incidence  and  observation.  This  procedure  is  used  by  Stevenson  [3] 
and  by  Asvestas  and  Kleinman  [4]  to  determine  the  low-frequency  electromagnetic  scattering 
from  perfectly  electrically  conducting  three-dimensional  bodies. 

However,  as  noted  in  [2]  this  definition  of  Rayleigh  scattering  cannot  be  used  for  two- 
dimensional  scattering  problems  since  the  scattered  field  in  general  does  not  have  a  conver¬ 
gent  expansion  in  powers  of  k.  The  two-dimensional  fields  cannot  be  expanded  in  powers 
of  k  because  the  two-dimensional  free-space  Green’s  function  is  a  Hankel  function  which 
has  a  branch  point  at  A:=0,  whereas  the  the  three-dimensional  free-space  Green’s  function  is 
analytic  in  k. 

In  the  literature  much  work  on  low-frequency  scattering  from  general  three-dimensional 
bodies  has  been  published  [2],  but  it  seems  that  only  van  Bladel  [5]  and  MacCamy  [6] 
have  done  work  on  low-frequency  scattering  from  general  two-dimensional  bodies.  In  [5]  the 
low-frequency  scattering  from  arbitrarily  shaped  dielectric  and  conducting  cylinders  with 
finite  cross  sections  is  considered.  For  transverse  magnetic  (TM)  polarization  [5]  uses  a 
time-harmonic  integral  equation  to  obtain  the  low-frequency  behavior  of  the  total  current 
from  which  the  low-frequency  scattered  field  is  calculated.  For  transverse  electric  (TE) 
polarization  an  expansion  of  the  current  in  integral  powers  of  k  is  assumed.  As  seen  from 
the  previous  discussion  and  from  the  exact  eigenfunction  solutions  for  the  circular  cylinder  [7, 
ch.2]  and  strip  [7,  ch.4]  such  an  expansion  does  not  exist.  However,  only  the  first  two  terms 
in  this  power  series  expansion  are  used  in  [5]  and  in  this  paper  it  will  be  shown  that  the  final 
TE  results  in  [5]  for  perfect  conductors  are  also  correct.  Reference  [6]  contains  the  formal 
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low-frequency  expansion  of  the  scattered  fields  from  two-dimensional  perfect  conductors  with 
finite  cross  section  in  cases  where  the  incident  field  can  be  expanded  in  even  powers  of  k. 
However,  no  explicit  expressions  for  the  expansion  coefficients  are  given. 

There  are  several  reasons  for  being  interested  in  two-dimensional  scattering  solutions. 
They  can  be  used  directly  as  approximate  solutions  to  certain  three-dimensional  scattering 
problems.  For  example,  the  field  close  to  the  middle  of  a  long  finite  rod  can  often  be 
well  approximated  by  the  field  from  the  corresponding  infinite  rod.  Also,  two-dimensional 
scattering  solutions  can  be  helpful  for  validating  computer  codes. 

Another  important  reason  for  being  interested  in  two-dimensional  solutions  is  that  they 
can  be  used  to  determine  three-dimensional  incremental  length  diffraction  coefficients.  These 
incremental  length  diffraction  coefficients  can  in  turn  be  used  to  determine  scattering  contri¬ 
butions  from,  for  example,  curved  ridges  (bumps)  and  channels  (dents)  that  have  constant 
cross  sections.  The  three-dimensional  incremental  diffraction  coefficients  can  be  found  di¬ 
rectly  from  the  corresponding  two-dimensional  far  fields  using  a  direct  substitution  approach. 
No  integration,  differentiation,  or  specific  knowledge  of  the  current  on  the  conductor  is 
needed.  This  direct  substitution  procedure  was  first  developed  by  Shore  and  Yaghjian  [8]  for 
planar  surfaces  and  then  extended  by  Hansen  and  Yaghjian  [9]  to  general  two-dimensional 
scatterers.  To  apply  this  procedure  one  must,  in  general,  be  able  to  evaluate  the  two- 
dimensional  far  fields  for  both  real  and  complex  angles  of  observation.  Therefore,  one  needs 
analytical  expressions  for  the  two-dimensional  far  fields. 

Analytical  expressions  for  two-dimensional  far  fields  from  single  cylindrical  bumps  and 
dents  can  also  be  used  in  conjunction  with  the  work  of  Twersky  [10]  to  calculate  the  scattering 
from  a  random  distribution  of  bumps  or  dents. 

The  purpose  of  this  report  is  to  determine  the  leading  terms  in  the  low-frequency  ex¬ 
pansions  of  the  scattered  electromagnetic  far  fields  for  the  following  three  types  of  two- 
dimensional  perfectly  electrically  conducting  scatterers: 

1.  An  arbitrarily  shaped  cylinder  with  finite  cross  section 

2.  An  arbitrarily  shaped  cylindrical  bump  (protuberance)  on  a  ground  plane 

3.  An  arbitrarily  shaped  cylindrical  dent  (indentation)  in  a  ground  plane. 

The  results  of  this  paper  apply  to  bumps  and  dents  that  are  continuously  lined  by  a 
conductor  and  thus  they  are  not  generally  valid  for  bumps  and  dents  that  contain  slits.  The 
scatterers  are  illuminated  by  plane  waves  propagating  in  directions  normal  to  the  axis  of 
the  cylinders,  and  the  two  polarization  cases,  TM  and  TE,  are  treated  separately.  Since 
the  scatterers  are  perfectly  conducting  the  low-frequency  solutions  for  normal  incidence  can 
be  immediately  generalized  to  obtain  low-frequency  solutions  for  oblique  incidence  (see,  for 
example,  [11,  sec.8.15]  or  [7,  ch.lj). 

As  mentioned  above,  the  unknown  field  for  two-dimensional  scatterers,  unlike  the  field 
for  three-dimensional  scatterers,  cannot  be  expanded  in  powers  of  k.  Instead  we  take  the 
following  approach:  The  Green’s  functions  and  incident  field  in  the  time-harmonic  integral 
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equation  for  the  current  is  expanded  for  small  k  to  obtain  an  integral  equation  for  the  low- 
frequency  current.  From  this  low-frequency  integral  equation  we  find  the  leading  terms  in 
the  low-frequency  expansion  of  the  current  and  calculate  the  low-frequency  far  field. 

The  low-frequency  expressions  for  the  scattered  fields  from  the  cylinders  with  finite  cross 
section  are  derived  in  Chapter  2.  For  the  cylindrical  bump  and  dent  the  scattered  fields 
are  written  as  the  sum  of  a  known  reflected  field  and  a  diffracted  field.  The  low-frequency 
expressions  for  the  diffracted  field  from  the  bump  and  dent  are  derived  in  Chapter  3  and 
Chapter  4,  respectively.  In  Chapter  5  the  low-frequency  expressions  are  verified  in  cases 
where  exact  time-harmonic  eigenfunction  solutions  exist  and  the  constants  are  calculated  for 
a  number  of  geometries.  A  summary  of  this  work  was  first  presented  in  [12]. 
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Chapter  2 

Cylinder  with  finite  cross  section 


The  cylinder  with  finite  cross  section  illuminated  by  the  plane- wave  field  is  situated 

in  a  rectangular  coordinate  system  shown  in  Figure  2.1.  The  scatterer  extends  uniformly  to 
infinity  in  the  -\-z  and  —z  directions  and  the  curve  that  describes  the  finite  cross  section  of 
the  scatterer  in  the  x  —  y  plane  is  denoted  by  S.  The  outward  normal  to  the  scatterer  is  n 
and  the  tangent  unit  vector  to5ist  =  ixn.  A  characteristic  dimension  of  the  scatterer  in 
the  X  —  y  plane  is  denoted  by  d.  In  addition  to  the  rectangular  coordinates  {x,y,z),  circular 
cylindrical  coordinates  (r,  <j>,  z)  given  by  x  =  r  cos  (f>,  y  =  r  sin  <f>,  and  z  =  z  will  also  be  used. 

The  incident  plane- wave  field  propagates  in  a  direction,  designated  by  4>\  normal  to  the 
2-axis.  As  usual,  the  scattered  field  is  defined  as  the  total  field  {E,H)  minus  the 

incident  field  {E\H'). 

Throughout  the  report  e-'"*  time  dependence  is  suppressed  in  all  the  time-harmonic  equa¬ 
tions. 


2.1  Transverse  magnetic  (TM)  polarization 

The  incident  electric  field  is  given  by 

E\f)  =  (2.1) 

where  k  is  the  propagation  constant.  The  low-frequency  current  is  determined  in  Section 

2.1.1  and  in  Section  2.1.2  it  is  integrated  to  get  the  low-frequency  far  field. 

2.1.1  Low-frequency  current  (TM) 

To  determine  the  low-frequency  far  field  we  first  determine  the  low-frequency  current.  As 
we  shall  see,  for  the  cylinder  with  finite  cross  section  for  TM  polarization  it  is  sufficient  to 
know  the  total  current  flowing  in  the  2— direction  on  the  cylinder.  To  find  this  current  we 
use  the  vector  potential  for  the  scattered  field 

A{f)  =  f  G{f,f')K{f')ds'  (2.2) 

Js 
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Figure  2.1:  Cylinder  with  finite  cross  section. 


from  which  the  scattered  fields  are  found  as 

xV  X  A,  H^  =  VxA. 

Here  G(f,  f  )  is  the  two-dimensional  free-space  Green’s  function 


G{f,f')  =  I  f  -  f' !)  (2.4) 

and  K{f)  =  zK^r)  is  the  current  on  the  cylinder.  Because  the  vector  potential  has  only  a 
^-component  and  the  free-space  Green’s  function  satisfies  the  scalar  homogeneous  Helmholtz 
equation  for  r  ^  f'  we  find  from  (2.2)  and  (2.4)  that 

E\if)  =  G(f,  f')K.(f')ds'  (2.5) 

which  is  a  relation  between  the  scattered  electric  field  and  the  current  on  the  cylinder. 

Equation  (2.5)  and  the  boundary  condition  of  zero  tangential  electric  field  on  the  con¬ 
ductor  yield  the  integral  equation  [13], [14] 


Em=}k^£G(f.f)K,{f')d,\  f€S  (2.6) 

where  the  bar  on  the  integral  sign  indicates  that  the  singularity  at  f  =  f'  is  excluded. 

From  the  small  argument  approximation  for  the  Hankel  function  we  find  that  for  f  ^  f’ 

G{f,r')  = -^lakd-\-0{{kd)°) ,  kd -k  0  (2.7) 
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where  d  is  a.  characteristic  dimension  of  the  cylinder. 

Because  the  incident  electric  field  in  (2.1)  for  low  frequencies  is  approximately  equal  to  i 
on  the  circumference  S  of  the  scatterer  it  is  assumed  that  the  total  current  for  low  frequencies 
is  non-zero,  that  is, 

/  A%(r')ds'  0.  (2.8) 

Js 

This  assumption  is  confirmed  by  the  eigenfunction  solutions  for  the  circular  cylinder  [7,  ch.2] 
and  strip  [7,  ch.4].  f 

Inserting  the  expansion  (2.7)  for  the  Green’s  function  into  the  integral  equation  (2.6)  and 
using  (2.8)  yields  the  low-frequency  behavior  of  the  total  current 

-  In  /  KAf')ds^  ~  1  ,  H  0.  (2.9) 

27r  V  t  Js 


2.1.2  Low-frequency  scattered  far  field  (TM) 

We  will  now  find  the  low-frequency  far  field  by  integrating  the  low-frequency  current  given 
in  (2.9).  From  the  asymptotic  expansion  of  the  Hankel  function  it  is  found  that 

G'(f,f')  ~  =[1 -I- jfc(i'cos</> -f  y'sin^)  4- Of(A:d)^)]  ,  r +oo.  (2.10) 

2v27r  vfcr 

W’hen  we  insert  (2.10)  into  the  integral  (2.5)  for  the  scattered  field  we  find  that  the  low- 
frequency  far  field  is  given  by 

Lp-j3ir/4  f 

Substituting  the  integral  of  the  current  from  (2.9)  into  the  field  expression  (2.11)  we  find  the 
final  low-frequency  expression  for  the  far  field  [5] 


E\{f)  ~ 


V  2  In  kd  -f/kr 


(2.12) 


Note  that  In  kd  ~  In  Akd  as  fed  0  when  A  is  a  positive  constant  so  that  d  can  be  any 
characteristic  dimension  of  the  cross  section  of  the  scatterer.  Furthermore,  it  is  noted  that 
the  low-frequency  scattered  field  does  not  depend  on  the  shape  of  the  finite  cross  section  but 
only  on  its  characteristic  dimension. 

The  expression  (2.12)  checks  with  the  low-frequency  results  for  the  circular  cylinder  [7, 
ch.2]  and  strip  [7,  ch.4]  obtained  from  exact  eigenfunction  solutions. 

It  should  be  noted  that  the  expression  in  (2.12)  for  the  far  field  cannot  be  obtained  from 
the  solution  to  the  purely  electrostatic  or  magnetostatic  problem  in  which  the  scatterer  is 
situated  in  the  static  field  E'\ic-o  or  H'\k=o,  respectively; 

Since  the  total  electrostatic  field  has  only  a  z— component  one  finds  from  Maxwell’s 
equations  that  it  must  be  a  constant  and  therefore  zero  because  it  is  zero  on  the  surface 
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of  the  scatterer.  Consequently,  the  scattered  electrostatic  field  equals  the  negative  of  the 
incident  field.  The  normal  component  of  this  electrostatic  solution  is  zero  everywhere  so 
that  the  charge  on  the  conductor  is  zero.  Moreover,  this  electrostatic  solution  gives  no 
information  about  the  current  on  the  conductor  and  thus  no  information  about  the  scattered 
time-harmonic  fields. 

As  part  of  the  solution  to  the  magnetostatic  equations  ( V  x  ^  =  0,  V  •  AT  =  0,  and 
h  •  //  =  0  on  5)  we  can  have  a  total  current  on  the  conductor.  However,  the  value  of  this 
current  is  not  determined  by  these  magnetostatic  equations. 


2.2  Transverse  electric  (TE)  polarization 

In  this  section  we  will  determine  the  low-frequency  far  field  in  the  case  of  TE  polarization. 
In  Section  2.2.1  we  determine  a  low-frequency  expansion  of  the  current  and  in  Section  2.2.2 
ve  integrate  this  low-frequency  expansion  to  get  the  far  field. 

The  incident  electric  and  magnetic  fields  are  here  given  by 

H'{f)  =  (2.13) 

and 

E'{f)  =  y^(sin<^'x  (2.14) 

2.2.1  Low-frequency  current  (TE) 

To  obtain  the  low-frequency  current  we  will  insert  the  low-frequency  expansions  for  the 
Green’s  function  and  the  incident  field  into  the  two-dimensional  magnetic  field  integral 
equation.  We  will  start  by  deriving  this  integral  equation.  Using  the  relations  (2.2)and 
(2.3)  invoVi  ■  the  vector  potential  and  the  fact  that 

V  X  (G(f,f')A-(f')|  =  (2.15) 

where  t  and  refer  to  the  tangent  -uid  normal  unit  vectors  in  Figure  2.1,  we  obtain 

(2.16) 

which  is  a  relation  between  the  current  and  the  scattered  magnetic  field. 

We  now  let  f  in  (2.16)  approach  a  point  on  S  and  use  the  result  from  [14,  App.j  of  an 
integration  near  the  singularity  at  f'  =  f  and  the  boundary  condition  =  —Kt  on  S  to  get 
the  magnetic  field  integral  equation  [14] 

1 

=  res.  (2.17) 
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The  bar  on  the  integral  sign  indicates  that  this  is  a  Cauchy  principal  value  integration.  To 
determine  the  low-frequency  expansion  of  the  current  we  use  the  small  argument  expansion 
of  the  Hankel  function  to  show  that 

=  g;;;G“(r,f')  +  0((W)"InW)  (2.18) 

where  G°(r,  r')  is  the  static  two-dimensional  free-space  Green’s  function  given  by 

=  (2.19) 

The  term  of  order  (kd)^  In  kd  in  (2.18)  is  non-singular  at  f  =  f'  ;  in  fact,  it  tends  to  zero  as 
r  — ►  r  . 

Letting  kd  0  in  the  integral  equation  (2.17)  and  using  the  expansion  (2.18)  we  find 
that  the  first  term  Kf  in  the  low-frequency  expansion  of  the  current  Kt  satisfies  the  integral 
equation 

1  =  £  -  iA-?(f),  1=  €  S.  (2.20) 

The  uniqueness  of  the  solution  to  this  integral  equation  has  been  proven  in  [15]. 

We  will  now  show  that  the  solution  to  the  integral  equation  (2.20)  is  related  to  a  magne¬ 
tostatic  problem:  A  static  magnetic  field  and  current  satisfy  the  relation  (2.16)  with  G{f,f') 
replaced  by  G®(f,  f').  Furthermore,  because  the  second  term  in  the  expansion  (2.18)  of  the 
derivative  of  the  Green’s  function  is  non-singular  at  f '  =  f  it  follows  that  the  integral  equa¬ 
tion  (2.17)  holds  in  the  static  case.  Equation  (2.20)  is  thus  the  integral  equation  for  the 
static  current  when  the  scatterer  is  situated  in  the  impressed  magnetostatic  field 

=  =  (2.21) 

Let  be  the  solution  to  this  static  problem.  Then  since  V  x  =  0  we  find  that  is 
a  constant.  Furthermore,  the  impressed  magnetostatic  field  in  (2.21)  satisfies  the  boundary 
condition  n  ■  =  0  on  the  conductor,  and  since  the  scattered  field  must  vanish  far  away 

from  the  scatterer  it  follows  that 


(2.22) 

We  can  therefore  write  the  current  A((r)  as 

/G(f)  =  -l-bA:[(f)  (2.23) 

where  KHf)  — +  0  as  — >  0.  We  have  hereby  obtained  the  first  term  in  the  low-frequency 

expansion  of  the  current  and  we  will  now  find  the  second  term. 

If  we  use  the  expansion  (2.23)  of  the  current,  the  expansion  (2.18)  of  the  Green’s  function, 
and  (2.20)  in  the  integral  equation  (2.17)  we  find  that  the  second  term  K}  in  the  low- 
frequency  expansion  of  the  current  Kt  satisfies  the  integral  equation 

jk{x  cos  <!>'  -\-y  sin  <f>')  =  f')ds'  -  f  £  S. 
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(2.24) 


Consequently,  the  second  term  is  of  order  kd  and  we  can  expand  the  current  as 

Kt{f)  =  -1  +  kdK^\r)  +  Kr{f)  (2.25) 

where  =  -qK}  and  >  0  as  fed  — »  0 

We  will  now  show  that  the  current  which  is  independent  of  kd,  can  be  found  from 
the  solution  to  the  electrostatic  problem  in  which  the  scatterer  is  situated  in  the  impressed 
electrostatic  field 

^'static  =  ^‘|fc=o  =  sin  <l>'  -y  cos  (2.26) 

Let  the  electrostatic  solution  be  denoted  by  E^.  Then  V  •  .E®  =  0,  that  is,  the  electrostatic 
field  is  solenoidal.  In  Appendix  A  it  is  shown  that  because  of  this  and  because  the  total 
charge  on  the  conductor  is  zero  the  field  ^EyX  4-  E^y  is  conservative.  We  can  therefore 
introduce  a  scalar  potential  so  that  —E^x  +  E^y  =  Defining  =  E^z  one 

finds  that 

^0  =  X  (2.27) 

and  the  electrostatic  field  is  thereby  written  in  terms  of  the  vector  potential  F°.  Because 
the  electrostatic  field  is  irrotational,  that  is,  V  x  =  0  it  follows  that  satisfies  Laplace’s 
equation  =  0.  Since  the  tangential  component  of  the  electric  field  is  zero  on  the  con¬ 

ductor,  satisfies  the  Neumann  boundary  condition  on  the  conductor.  Hence,  satisfies 
the  same  differential  equation  (Laplace’s  equation)  and  the  same  boundary  conditions  on  the 
conductor  (Neumann  type)  as  the  z-component  of  a  magnetostatic  field.  If  we  furthermore 
require  that  the  potential  Ff*  for  the  scattered  electrostatic  field  tend  to  zero  at  infinity  we 
find  that  F^  satisfies  the  integral  equation  (2.17)  with  Hl,G,  and  Kt  replciced  by  F^,G°, 
and  —  F®,  respectively.  The  impressed  electrostatic  field  given  in  (2.26)  corresponds  to  the 
impressed  vector  potential 


i^(r)  =  z  cos  -1- y  sin  -H  G  (2.28) 

where  G  is  a  constant.  The  integral  equation  for  F^  is  thus 

z cos <l>'  +  ysm<l>'  +  C  =  -£  Ef(f')-^G°(f,f')ds'  +  ^i^(f),  f  €  5.  (2.29) 

Using  the  result  from  [14,  App.J  of  the  integration  near  the  singularity  at  f'  =  r  and  the 
divergence  theorem  in  the  region  bounded  by  S  one  finds  that 

C  =  -£^G”(f,f')C*'+lc,  res.  (2.30) 

Consequently,  the  constant  function  C  on  the  left  hand  side  of  the  integral  equation  (2.29) 
only  adds  the  constant  G  to  the  solution  to  (2.29)  with  G  =  0.  Because  we  may  add  a 
constant  to  the  potential  F°  without  changing  the  field  E°  we  can  let  G  =  0  in  the  integral 
equation  (2.29). 
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By  comparing  the  integral  equations  (2.29)  with  C  =  0  and  (2.24)  and  noting  that  the 
solutions  are  unique  one  finds  Kt^(r)  =  —  j^(r),  r  £  S,  which  completes  the  proof  that  the 
second  term  in  the  low-frequency  expansion  of  the  current  can  be  found  from  the  electrostatic 
solution.  The  low-frequency  expansion  of  the  current  is  thus 

K,(f)  =  -l-jkF;(f)  +  K"(f)  (2.31) 

where  satisfies  the  electrostatic  integral  equation  (2.29)  and  0  bls  kd  —*  0. 

In  this  Section  we  have  shown  that  the  first  and  second  term  in  the  TE  low-frequency 
expansion  of  the  current  are  of  order  (fcd)°  and  fed,  respectively.  Since  the  second  term  in 
the  low-frequency  expansion  (2.18)  of  the  derivative  of  the  free-space  Green’s  function  is  of 
order  {kd)^  In  kd  it  is  seen  from  the  integral  equation  (2.17)  that  the  third  term  in  the  low- 
frequency  expansion  of  the  current,  in  general,  will  be  a  function  of  kd  and  In  kd.  Therefore, 
as  noted  in  the  Introduction,  the  current  cannot  be  expanded  in  a  power  series  in  kd. 

2.2.2  Low-frequency  scattered  far  field  (TE) 

We  will  now  calculate  the  far  field  scattered  by  the  current  given  in  the  expansion  (2.31). 
Using  the  asymptotic  expansion  of  the  Hankel  function  one  finds  that 

~G{f,f')  ~  k^-j=^-^[l  +  jk{x' cos  <f>  +  y' sin  <f>)  +  0{{kdy)]h'  ‘f  (2.32) 

on'  2v2?r  Vkr 

as  r  ^  -foo.  Inserting  (2.32)  and  the  current  expansion  (2.31)  into  the  expression  (2.16)  for 
the  scattered  magnetic  field  one  finds  that  the  scattered  field  to  order  (kd)^  is  given  by 


p-jjr/4  -jkr  i  .  ^ 

H>(f)  ^  _(A:d)2l— 1— — [r  •  /  {x'cos4>  +  y'sm<f>)h'ds'  +  f-  /  F^{f')h'ds'].  (2.33) 

2v27r  y/kr  d^  Js  JS 


Applying  the  divergence  theorem  one  finds  that  the  first  term  in  the  brackets  of  (2.33)  equals 
the  area  bounded  by  S,  that  is,  the  area  of  the  cross  section  of  the  scatterer.  Integrating  the 
second  term  in  the  brackets  of  (2.33)  by  parts,  aad  using  that 


C7°(r-) 


=  n-B 


r  £  S. 


(2.34) 


where  cr°  is  the  electrostatic  charge  on  the  conductor,  we  find  that  the  scattered  field  is  given 

by 


p-J)r/4  p-jkr  i  rr  . 


S  c 


(2.35) 


where  /I5  is  the  area  of  the  cross  section  of  the  scatterer.  From  this  expression  it  is  seen 
that  the  scattered  field  consists  of  a  magnetic  dipole  in  the  z— direction  (the  term  with  As) 
plus  an  electric  dipole  in  the  x  —  y  plane  (the  integral  term)  [11,  Sec.3.8]. 

This  expression  for  the  scattered  field  may  be  rewritten  in  a  form  that  is  more  convenient 
for  the  following  section  which  deals  with  the  cylindrical  bump.  From  the  expression  (2.26) 
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for  the  impressed  electrostatic  field  it  follows  that  the  total  electrostatic  charge  <t“  can 

be  written  as 

sin  cos  <f>*  (2.36) 


where  and  <r°®'  are  the  electrostatic  charges  in  the  electrostatic  problems  where  the 
scatterer  is  situated  in  the  impressed  electrostatic  fields  x  and  y,  respectively.  Using  the 
expression  (2.36)  for  cr®  and  the  dipole  moment  reciprocity  theorem 


/  =  /  a^(r^)x'd3' 

J  S  J  s 


proven  in  Appendix  B  the  scattered  magnetic  far  field  may  be  written  as 


+C2  cos  cos  <j)'  +  C3  sin((^  + 


(2.37) 


(2.38) 


where 


<T®*(f) 


f  f 

xds,  €2=  - -yds,  Ca  =  -  / 

Js  t  Js 


<7®*'(r) 


(2.39) 


We  will  now  present  integral  equations  that  determine  the  electrostatic  charges  <t®*  and 
An  integral  equation  for  can  be  derived  from  the  scalar  potential  0  for  the  electro¬ 
static  field.  The  potential  ^  satisfies  Laplace’s  equation  and  the  potential  for  the  impressed 
electrostatic  field  is  given  by  i/>'  =  -i  -hC,  where  C  is  a  constant.  Noting  that  the  potential 
for  the  scattered  field  is  given  by  V’*  =  !$  and  that  the  total  potential  is  constant 

on  S,  one  obtains  an  integral  equation  involving  an  unknown  constant.  This  constant  is  eval¬ 
uated  at  an  arbitrary  observation  point  fo  =  (xo,  J/o)  €  5  so  that  a®*  satisfies  the  integral 
equation 

/  [G"(f,f')  -  =  I  _  ar„,  f^S.  (2.40) 

J  s  c 

In  [15]  it  is  shown  that  the  integral  equation  (2.40)  along  with  the  condition  that  the  total 
charge  on  the  conductor  is  zero,  that  is,  <r®*<Zs  =  0,  determine  the  electrostatic  charge  cr®’^ 
uniquely.  Similarly,  cr®*'  is  determined  by 


”('■0. Wo,  reS. 


(2.41) 


along  with  the  condition  that  the  total  charge  on  the  conductor  is  zero,  that  is,  jg  a^ds  =  0. 

Alternatively,  the  three  constants  Ci,C2,  and  C3  in  (2.39)  can  be  written  in  terms  of 
vector  potentials  for  the  electrostatic  fields.  From  the  integral  equation  (2.29)  with  6'  =  0  it 
is  seen  that  we  may  introduce  vector  potentials  F®*  and  satisfying 


(2.42) 
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and 


y  =  -£  ^G“(f , r')i^(r-')*'  +  1  J^(f ),  f  €  5 


(2.43) 


SO  that 

F®  =  cos  +  sin  (2.44) 

The  integral  equations  (2.42)  and  (2.43)  are  of  the  same  form  as  the  time-harmonic  two- 
dimensional  magnetic  field  integral  equation  and  they  should  therefore  be  easy  to  solve 
numerically  when  no  part  of  the  scatterer  is  flat  [16,  p.l68].  (For  flat  scatterers,  (2.40)  and 
(2.41)  can  be  used.)  Using  the  relation  (2.34)  between  <7°  and  and  the  relations  (2.36) 
and  (2.44)  one  finds  that 

^  =  res  (2.45) 

and 

=  res.  (2.46) 

Inserting  the  relations  (2.45)  and  (2.46)  into  the  expressions  (2.39)  for  the  constants  Ci,C2, 
and  C3  and  integrating  by  parts  one  finds 

which  gives  the  three  constants  in  terms  of  the  vector  potentials  F®*  and  F^. 

We  have  now  given  the  exact  analytical  expressions  for  the  first  terms  in  the  TM  and 
TE  low-frequency  scattered  far  fields  from  the  cylinder  with  finite  cross  section.  We  found 
that  these  expressions  were  completely  determined  by  Ccilculating  three  constants  that  only 
depended  on  the  shape  of  the  cross  section  of  the  cylinder.  These  three  constants  are  in  turn 
found  from  the  electrostatic  solutions  for  x  and  y  directed  impressed  electrostatic  fields. 
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Chapter  3 

Cylindrical  bump  on  an  infinite 
ground  plane 


In  this  chapter  the  low-frequency  expressions  are  derived  for  the  field  diffracted  by  a  cylin¬ 
drical  bump  on  a  ground  plane  illuminated  by  a  plane  wave.  Both  the  bump  and  the  ground 
plane  are  perfectly  conducting.  The  bump  extends  uniformly  to  infinity  in  the  +2  and  —2 
directions  and  the  curve  that  describes  the  cross  section  of  the  bump  in  the  x  —  y  plane 
is  denoted  by  B  as  shown  in  Figure  3.1.  The  outward  normal  to  the  bump  is  n  and  the 
tangent  unit  vector  to  B  is  t  =  2  x  n.  The  equation  for  the  ground  plane  is  y  =  0  and  the 
intersections  between  the  bump  and  the  ground  plane  are  given  by  y  =  0,  a:  =  ±d. 

The  incident  field  {E\H')  and  the  cylindrical  coordinates  {r,<j>,z)  are  the  same  as  in 
Section  2.  The  field  reflected  in  the  ground  plane  y  =  0  when  there  is  no  bump  is  denoted 
by  (£’■,//’’■).  We  write  the  total  field  as  {E,H)  =  {E\H')  +  -|-  where 

is  by  definition  the  diffracted  field. 

It  will  be  shown  that  the  diffracted  field  is  the  scattered  field  in  an  equivalent  scattering 
problem  where  the  scatterer  is  the  bump  plus  its  image  in  the  ground  plane,  and  the  incident 
field  is  {E'  +  E^,H'  -}-  H^)  as  shown  in  Figure  3.2.  The  curve  that  describes  the  image  of  the 
bump  in  the  y  =  0  plane  is  denoted  by  B,-  and  the  cross  section  of  the  equivalent  scatterer 
is  therefore  enclosed  by  5  U  Bi. 

By  considering  the  TM  and  TE  polarizations  separately  it  is  eaisily  seen  that  the  scattered 
field  in  the  equivalent  scattering  problem  satisfies  the  boundary  conditions  on  the  ground 
plane.  Since  {E'  -f  H'  -j-  H'")  also  satisfies  these  boundary  conditions  it  is  clear  from  the 
uniqueness  of  the  solution  that  the  scattered  field  in  the  equivalent  scattering  problem  is  the 
diffracted  field  in  the  bump  scattering  problem.  This  method  of  images  for  constructing  an 
equivalent  scattering  problem  was  first  used  by  Rayleigh  [17]  when  he  solved  for  scattering 
from  the  semi-circular  bump. 
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Figure  3.1:  Cylindrical  bump  on  a  ground  plane. 
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3.1  TM  polarization 

We  first  consider  the  TM  case  where  the  incident  electric  field  is  given  by  (2.1)  and  the 
incident  magnetic  field  is  given  by 

H\f)  =  ^(-  sin  <f>'x  +  cos  ^iyyk(rccei,>+ysin4>')  (3 

The  reflected  electric  field  is  given  by  (2.1)  with  z  and  y  replaced  by  —z  and  — y,  respectively. 
Similarly  the  reflected  magnetic  field  is  given  by  (3.1)  with  y  and  y  replaced  by  — y  and  — y, 
respectively. 

If  we  use  the  formula  (2.12)  for  the  scattered  far  field  and  eidd  the  contributions  from  the 
incident  field  E'  and  the  reflected  field  we  get  zero  because  the  two  contributions  cancel. 
Consequently,  the  diffracted  field  is  zero  to  order  and  we  cannot  use  the  results  from 
Section  2.1  because  we  need  higher  order  terms.  We  will  therefore  start  by  investigating  the 
low-frequency  current  on  the  equivalent  scatterer  given  by  B  U  S,. 

3.1.1  Low-frequency  current  (TM) 

From  the  integral  equation  (2.6)  it  is  found  that  the  current  on  the  equivalent  scatterer 
satisfies  the  integral  equation 

+  £;(f )  =  /  G(r,  re  BUB,.  (3.2) 

From  the  small  argument  approximation  of  the  Hankel  function  it  is  found  that 

G(f ,  f')  =  G«(r,  r-)  -  i  In  ^  i  -  ^  +  0((  In  kd)  (3.3) 

where  G°(f,f')  is  the  static  free-space  Green’s  function  in  (2.19),  7  is  Euler’s  constant,  and 
the  last  term  is  non-singular  at  f  =  f',  in  fact  it  tends  to  zero  as  f  — »  f'.  Because  El  -H  El 
is  an  odd  function  of  y  and  the  equivalent  scatterer  is  symmetric  around  the  y  =  0  plane  it 
follows  that  the  current  is  also  an  odd  function  of  y,  that  is, 

Kz{x,  y)  =  -y).  (3.4) 

Consequently,  the  total  current  flowing  along  the  equivalent  scatterer  is  zero,  that  is, 

f  I<,{f)ds  =  0.  (3.5) 

JBuBi 

If  we  expand  El  +  El  in  powers  of  k,  and  use  the  expansion  (3.3)  of  the  Green’s  function  along 
with  (3.5)  in  the  integral  equation  (3.2),  we  find  that  the  first  term  in  the  low-frequency 
expansion  of  the  current  satisfies  the  integral  equation 

2ysin<^‘  =  y^/  G^{f,f')K°{f')ds\  feBuB..  (3.6) 
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The  uniqueness  of  the  solution  to  this  integral  equation  has  been  proven  in  [15]. 

We  will  now  show  that  K°  is  the  current  in  the  magnetostatic  scattering  problem  with 
the  equivalent  scatterer  situated  in  the  impressed  magnetostatic  field 


SLtic  =  S'\k=o  + 


The  magnetostatic  current  in  this  scattering  problem  with  the  equivalent  scatterer  is  denoted 
by  and  a  vector  potential  for  the  scattered  magnetostatic  field  is  given  by  (2.2)  with  G 
replaced  by  and  K  replaced  by  zK°^.  From  the  boundary  condition,  n  •  =  0,  for  the 

magnetostatic  field  we  find  that  the  total  magnetostatic  vector  potential  must  be  constant 
on  the  conductor,  and  we  therefore  require  that  it  be  zero  there.  From  the  expression  (3.7) 
for  the  impressed  magnetostatic  field  it  is  seen  that  the  impressed  vector  potential  must 
be  — 2j/y^sin^‘  +  C  where  C  is  a  constant,  which  is  undetermined  at  this  point  of  the 
derivation.  An  integral  equation  for  is  thus 

2y . sin  -f  G°{f,  f')/r®  (f')ds',  f  €  B  U  B..  (3.8) 

y  JBuBi 

Because  the  impressed  magnetostatic  field  is  parallel  to  x  and  the  equivalent  scatterer  is 
symmetric  around  the  plane  y=Q  the  current  excited  by  this  impressed  field  must  be  an  odd 
function  of  y.  Therefore,  the  integral  in  (3.8)  must  be  an  odd  function  of  y  and  the  constant 
C  must  be  zero. 

If  we  compare  (3.8)  with  C  =  0  to  the  integral  equation  (3.6)  for  the  current  we  see 
that  A'fj  =  A'®  and  we  have  shown  that  A"°  is  the  current  in  the  magnetostatic  problem. 
The  leading  term  in  the  low-frequency  expansion  of  the  current  is  thus  of  order  (A:d)°.  From 
the  integral  equation  (3.2)  it  follows  that  the  next  term,  in  general,  is  a  function  of  kd  and 
In  kd  so  that  the  current  cannot  be  expanded  in  powers  of  kd. 

It  is  convenient  to  define  a  current  which  is  the  magnetostatic  current  when  the 
equivalent  scatterer  is  situated  in  the  impressed  magnetostatic  field  x.  From  the  expression 
(3.7)  for  the  impressed  magnetostatic  field  and  from  the  integral  equation  (3.6)  it  follows 
that  the  current  A'°®  is  determined  by  the  integral  equation 


-y=i  rgBUBj 

JBuBi 


and  the  condition  (3.4)  that  it  must  be  an  odd  function  of  y.  We  then  can  write  the  low- 
frequency  expansion  of  the  current  as 


Kz  =  —24/  — sin  <^' 
VM 


Kf  +  k: 


(3.10) 


where  A'”®  is  an  odd  function  of  y  and  satisfies  the  magnetostatic  integral  equation  (3.9), 


and  Kl  — >  0  as  — ♦  0. 
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3.1.2  Low-frequency  diffracted  far  field  (TM) 

We  will  now  calculate  the  low-frequency  diffracted  far  field  from  the  bump,  which  is  the 
scattered  far  field  in  the  equivalent  scattering  problem  of  the  bump  and  its  image. 

Using  the  electric  field  expression  (2.5)  and  the  far-held  approximation  (2.10)  to  the 
Green’s  function  we  find  that  the  diffracted  far  field  is 

Et{r)  ~  -jk^  -  I  [1  +  jkix' cos  +  y' sin  <f>)  +  0{{kdf)]K^{f ')ds' .  (3.11) 

2v2ir  vkr  V  c  JbuB, 

Because  the  current  is  an  odd  function  of  y  and  because  the  equivalent  scatterer  is  symmetric 
around  the  ground  plane  y  =  0  we  find  that 

/  xK^{f)ds  =  0,  /  yKz{r)ds  =  2  /  yK^{f)ds.  (3.12) 

JBuBi  JBuBi  Jb 

Furthermore,  using  the  fact  that  the  total  current  is  zero  (3.5)  together  with  the  expansion 
(3.10)  of  the  current  we  find  that  the  diffracted  far  field  (3.11)  to  order  (kd)^  is  given  by 

E'^{f)  ~  —{kd)^e~^^^‘*J - Bosin<^  sin<^'  (3.13) 

V  Vfcr  « 

where 

Bo=  f  yKf{f)ds.  (3.14) 

J  B 

The  magnetostatic  current  AT”®  is  an  odd  function  of  y  and  satisfies  the  static  integral 
equation  (3.9). 

From  the  formula  (3.13)  it  is  seen  that  the  diffracted  field  is  the  field  of  a  magnetic  dipole 
in  the  x-direction  [11,  sec. 3. 8]. 

3.2  TE  polarization 

In  the  TE  case  the  incident  magnetic  and  electric  fields  are  given  in  (2.13)  and  (2.14). 
Because  the  contributions  from  the  incident  and  reflected  fields  in  the  TE  case  do  not  cancel 
we  can  use  the  results  from  Section  2.2  directly  so  we  do  not  have  to  investigate  the  current 
on  the  equivalent  scatterer. 

3.2.1  Low-frequency  diffracted  far  field  (TE) 

If  we  let  Ab  denote  the  area  of  the  region  bounded  by  B  and  the  line  y  =  0,  |x|  <  d  (see 
Figure  3.1)  then  the  diffracted  far  field  is  given  by 

Hi(f)  =  (3.15) 
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where  //*  is  the  scattered  field  given  in  (2.38).  The  first  term  in  (3.15)  is  the  contribution  from 
{E',H')  and  the  second  term  is  the  contribution  from  Because  of  the  symmetry 

of  the  scatterer  we  find  that 


y<T^{r)ds  =  2  /  ya^{r)ds 

J  B 


(3.16) 


where  and  are  the  electrostatic  charges  in  the  electrostatic  problems  with  the  equiv¬ 
alent  scatterer  situated  in  the  impressed  electrostatic  fields  x  and  y,  respectively. 

With  substitution  of  (3.16)  into  the  field  expression  (2.38),  equation  (3.15)  becomes 


Hi(f)  ~  -(W)=e-"'‘,/i^l|/le  +  cos,Acos,^-  /  (3.17) 

V  ^  y/kr  d^  Jb  e 

The  diffracted  field  therefore  consists  of  the  contribution  from  a  magnetic  dipole  in  the 
2— direction  and  an  electric  dipole  in  the  y-direction  [11,  sec. 3. 8]. 

Alternatively,  using  (2.45)  the  integral  in  (3.17)  may  be  written  in  terms  of  the  vector 
potential  satisfying  (2.42)  with  S  replaced  by  B  U  Bi. 

W/e  will  now  show  that,  remarkably,  the  integral  in  the  expression  (3.17)  for  the  TE 
diffracted  field  equals  the  negative  of  the  constant  Bq  in  (3.14)  occurring  in  the  expression 
(3.13)  for  the  TM  diffracted  field.  To  do  this  we  introduce  a  scalar  potential  for  the  elec¬ 
trostatic  field  and  use  the  fact  that  the  impressed  scalar  potential  must  be  C  —  y  where  C 
is  a  constant.  Furthermore  the  boundary  condition  of  zero  tangential  electric  field  requires 
that  the  total  scalar  potential  be  constant  on  the  equivalent  scatterer,  and  we  find  that  the 
charge  therefore  must  satisfy  the  integral  equation 

y-C=-f  G^{f,f')^^^ds\  f€B\JBi  (3.18) 

JbuB,  € 

where  C  is  undetermined  at  this  point  of  the  derivation. 

Because  the  impressed  electrostatic  field  is  symmetric  around  the  ground  plane  y  =  0  the 
electrostatic  charge  must  be  an  odd  function  of  y.  Consequently,  the  integral  in  (3.18) 
must  be  an  odd  function  of  y  and  the  constant  C  must  be  zero. 

Comparing  the  integral  equation  (3.18)  with  C  =  0  to  the  integral  equation  (3.9)  for  the 
current  one  finds  that  and  therefore  the  TE  diffracted  field  is  given  by 


H^ir)  ~ 


1 

—^^[Ab  -  Bo  cos  4>  cos  (f>'] 
y/kr  a 


(3.19) 


where  Bq  is  given  in  (3.14)  or  (3.17). 

We  have  now  derived  the  exact  expressions  for  the  leading  terms  in  the  low-frequency 
expansions  of  the  diffracted  far  fields  from  arbitrarily  shaped  bumps.  Our  general  analysis 
has  proven  the  rather  remarkable  result  that  the  constant  in  the  TM  magnetic  dipole  term 
equals  the  negative  of  the  constant  in  the  TE  electric  dipole  term  for  arbitrarily  shaped 
bumps.  In  the  following  chapter  we  will  show  that  this  result  also  holds  for  diffraction  from 
a  two-dimensional  dent  in  a  ground  plane. 
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Chapter  4 

Cylindrical  dent  in  an  infinite  ground 
plane  ^ 


In  this  chapter  we  derive  the  low-frequency  expressions  for  the  fields  diffracted  by  a  cylindrical 
dent  in  a  ground  plane  illuminated  by  a  plane  wave.  Both  the  dent  and  the  ground  plane  Me 
perfectly  conducting.  The  dent  extends  uniformly  to  infinity  in  the  +z  and  —z  directions 
and  the  curve  that  describes  the  cross  section  of  the  dent  in  the  x  —  y  plane  is  denoted  by 
D  as  shown  in  Figure  4.1.  The  ground  plane  is  given  by  y  =  0  and  the  dent  intersects 
the  ground  plane  at  y  =  0,  x  =  ±d.  The  line  segment  given  by  y  =  0,  |x|  <  d  that  caps 
the  dent  is  called  the  aperture  and  is  denoted  by  A.  The  image  of  D  with  respect  to  the 
ground  plane  y  =  0  is  denoted  by  £),,  and  the  part  of  the  upper  half  plane  (y  >  0)  which  is 
outside  the  image  £),  is  called  F.  The  normal  to  D  is  n  and  the  tangent  unit  vector  to  D  is 
t  =  z  X  n.  Normal  and  tangent  unit  vectors  to  A  are  y  and  — x,  respectively.  The  incident 
field  the  reflected  field  the  diffracted  field  {E'^,H'^),  and  the  cylindrical 

coordinates  (r,  z)  are  defined  as  in  the  previous  chapter. 

Although  the  dent  problem  is  very  similar  to  the  bump  problem,  the  analysis  of  the  dent 
is  much  more  complicated  and  thus  warrants  a  separate  analysis.  The  main  reason  for  the 
complication  is  that  the  dent  problem,  unlike  the  bump  problem,  cannot  be  reduced  using 
image  theory  to  an  equivalent  problem  that  involves  only  a  cylinder  of  finite  cross  section 
without  a  ground  plane.  One  therefore  has  to  deal  with  integral  equations  more  complicated 
than  the  usual  two-dimensional  electric  and  magnetic  field  integral  equations.  It  is  essential 
to  the  derivation  that  the  integral  equation  involves  only  fields  and  currents  in  a  region 
whose  maximum  dimension  tends  to  zero  as  the  characteristic  dimension  of  the  dent  tends 
to  zero.  We  found  that  the  most  convenient  integral  equations  for  this  purpose  were  coupled 
integral  equations  first  derived  by  Asvestas  and  Kleinman  [18].  In  Hansen  and  Yaghjian  [19] 
these  integral  equations  have  been  derived  straightforwardly  by  means  of  the  Stratton-Chu 
formulas  [20].  Here  they  will  be  derived  as  in  [18]  by  means  of  Green’s  identities. 
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4.1 

Figure  4.1:  Cylindrical  dent  in  a  ground  plane. 

Transverse  magnetic  (TM)  polarization 

The  incident  electric  and  magnetic  fields  are  here  given  by 

£»^f)  = 

(4.1) 

and 

//•(f)  =  ./^(-sin<^‘i  -1-  cos^‘j))e-’*<^'“^‘+>'""^'). 

(4.2) 

The  reflected  electric  field  is  given  by  (4.1)  with  z  and  y  replaced  by  — z  and  —y,  respectively. 
Similarly,  th«-  reflected  magnetic  field  is  given  by  (4.2)  with  y  and  y  replaced  by  —y  and  —y, 
respectively. 

4.1.1  Integral  equation  for  the  current  (TM) 

We  will  start  by  giving  a  short  derivation  of  Asvestas  and  Kleinman’s  [18]  coupled  TM 
integral  equations  for  the  current  in  the  dent  and  the  magnetic  field  in  the  aperture.  We 
need  the  details  of  this  derivation  to  derive  our  low-frequency  results. 

Let  Si{R)  be  the  part  of  the  z— axis  given  by  |x|  <  R  with  R  >  d  euid  let  S2{R)  be  the 
part  of  the  circle  r  =  R  which  is  in  the  region  y  >  0.  Furthermore,  let  ns  be  that  normal 
unit  vector  to  St(R)  U  S2(R)  which  points  into  the  region  enclosed  by  Si(R)  U  S2(R)  shown 
in  Figure  4.2.  Let  Go  be  the  Dirichlet  Green’s  function  for  the  region  j/  >  0,  that  is, 

Goir,  f')  =  G(f,  r)  -  G(f,  f')  (4.3) 
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Figure  4.2:  Half-circle  in  the  upper  half- plane  y  >  0. 

where  G(r,  f')  is  the  two-dimensional  free-space  Green’s  function  in  (2.4)  and 

fi=:xx  —  yy  (4-4) 

is  the  image  in  the  ground  plane  of  the  point  f.  If  we  use  Green’s  second  identity  and  the 
fact  that  the  diffracted  field  satisfies  the  scalar  homogeneous  Helmholtz  equation  we  find 
that  the  diffracted  field  satisfies 

EHf)  =  -  /  lGD(r.  -  Ei(f)VGo(f,  r’))  •  h'^ds',  j,  >  0.  (4X) 

Because  of  the  boundary  conditions  for  Gq  and  on  the  ground  plane,  (4.5)  becomes 

EHr)  = 

+  y>a.  (4.6) 

Using  the  asymptotic  relations  for  Gp  and  that  E^  satisfies  the  two-dimensional  radiation 
condition  [11,  sec.1.34]  one  may  show  that  the  last  integral  in  (4.6)  tends  to  zero  as  R  oo. 
Thus  (4.6)  implies 

=  j^Ei(f)^GD{f,f)d3',  !/>0.  (4.7) 
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Using  the  property  of  the  Dirichlet’s  Green’s  function  that  ^Goir,  r')  =  2^G(r,  r'),  r'  E  A 
and  that  £’‘(r)  +  =  0,  f  e  A,  (4.7)  becomes 

E.(f)  =  2jf  £,(f')^G(f,f')*'+  f;;(f)  +  EHr),  y  >  0  (4.8) 

which  is  an  expression  for  the  total  field  outside  the  dent  in  terms  of  the  field  in  the  aperture. 
Taking  the  curl  of  (4.8)  to  get  the  magnetic  field,  we  find 

-  +  HKf)],  j>0.  (4.9) 

Because  the  integral  in  this  equation  is  an  even  function  of  y  we  have 

-  jk^[HUfi)  +  H:(fi)l  feinl(AUD)  (4.10) 

where  f,  is  given  in  (4.4)  and  int(A  U  D)  is  the  region  bounded  by  A\J  D.  To  get  rid  of 
the  integral  involving  the  second  derivative  of  the  Green’s  function  we  will  derive  another 
relation  involving  this  second  derivative  and  then  combine  these  two  relations. 

If  we  recall  that  Ez  satisfies  the  homogeneous  Helmholtz  equation  and  apply  Green’s 
second  identity  in  the  region  bounded  by  A  U  D  we  find  that 

Ez{r)=  !  [G{f,r)—Ez{r)-Ez{f')-^G{f,f')]ds\  f  G  mt(/l  U  T>).  (4.11) 

Ja\jd  on'  On' 

Taking  the  curl  of  (4.11)and  applying  the  boundary  condition  of  zero  tangential  electric  field 
in  the  dent  produces  the  following  equation  for  the  x— component  of  the  magnetic  field, 

+  f  Ez{f')-^-rG{f,f')ds',  f€int{AUD).  (4.12) 

Ja  Oy'* 

Insert  the  integral  in  (4.10)  involving  the  second  derivative  of  the  Green’s  function  into  (4.12) 
to  get 

-kfy  L 

+Hiifi)  +  H:{fi)],  feint{AuD).  (4.13) 
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We  now  let  f  approach  the  aperture  and  use  the  results  of  the  integration  near  the  singularity 
in  [14,  App.],  the  induction  law  of  Maxwell’s  equations,  the  following  expression  for  the 
current  on  D, 

and  the  relation 

Q 

^G(r,f')  =  0,  y  =  y'  =  0,  x  ^  x'  (4.15) 

to  get  from  (4.13) 

/d  f  e 

which  is  a  relation  between  the  current  in  the  dent  and  the  magnetic  field  in  the  aperture. 

To  determine  the  current  in  the  dent  we  need  one  more  relation  of  this  kind.  To  obtain 
this  we  apply  Green’s  second  identity  in  the  region  bounded  hy  A\J  D  and  find  that 

E,(f)  =  j^GN(f,f')-~E.(?)ds' 

t 

+  2  j^G{f,f')-^Ez{r')ds',  feint{AUD)  (4-17) 

where  Gi^  is  the  Neumann  Green’s  function  for  the  region  y  >  0,  that  is, 

G/v(f,f')  =  G(f,f')  +  G(f,f,-),  (4.18) 

Let  f  approach  D  in  (4.17),  and  use  the  boundary  condition  of  zero  tangential  electric  field 
in  the  dent  along  with  the  expression  (4.14)  for  the  current,  and  the  induction  law  to  get 
the  second  integral  equation 

/  Gyv(f,  f')K,{f')ds'  +  2  [  Gif,  f')H,if')ds'  =  0,  f  £  D.  (4.19) 

Jd  Ja 

The  relations  (4.16)  and  (4.19)  constitute  a  pair  of  coupled  integral  equations  that  determine 
the  current  Kz  in  D  and  the  magnetic  field  Hx  in  A. 

4.1.2  Low-frequency  current  (TM) 

We  will  now  use  the  coupled  integral  equations  (4.16)  and  (4.19)  to  determine  the  first  term 
in  the  low-frequency  expansion  of  the  current  in  the  dent. 

Using  the  low-frequency  expansion  (2.18)  of  the  derivative  of  the  Green’s  function  in  the 
integral  equation  (4.16)  we  find  that 

Jd  "  e  A  (4.20) 
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where  and  are  the  first  terms  in  the  low-frequency  expansions  of  the  current  in  the 
dent  and  the  magnetic  field  in  the  aperture,  respectively.  is  the  static  two-dimensional 
free-space  Green’s  function  given  in  (2.19).  By  using  the  low-frequency  expansion  of  the 
free-space  Green’s  function  (3.3)  and  the  corresponding  expansion  of  the  Neumann  Green’s 
function  in  the  integral  equation  (4.19)  we  find  that 


/  K°Ar')ds'  -I-  /  Hl{f')ds'  =  0 
Jd  Ja 


(4.21) 


and 


£  +  2  £  G«(r,f')/f.“(r')(is'  =  0, 


f  e  D. 


(4.22) 


The  relations  (4.20),  (4.21),  and  (4.22)  constitute  a  set  of  coupled  integral  equations  that 
determine  the  low-frequency  current  on  D.  In  [15]  it  is  shown  that  these  three  coupled 
integral  equations  have  a  unique  solution. 

We  will  now  show  that  these  integral  equations  are  the  integral  equations  for  the  mag¬ 
netostatic  current  when  the  scatterer  is  situated  in  the  impressed  magnetostatic  field 


SLuc  = 

We  consider  the  magnetostatic  problem  in  which  the  magnetostatic  solution  satisfies 
n  ■  Hi  =  0  and  h  x  =  —Ki  on  the  conductor.  Furthermore,  it  cein  be  proven  that  these 
boundary  conditions  make  the  total  static  field  zero  below  the  ground  plane  and  dent.  Since 
the  total  magnetostatic  field  is  zero  below  the  ground  plane  and  dent  it  is  convenient  to  divide 
the  total  field  into  three  components  as  Hi  =  Hltatic  +  ^Itatic  +  ^static  where  Hhatic 
scattered  field  when  there  is  no  dent  and  H'^tatic  definition  the  diffracted  magnetostatic 
field. 

We  introduce  a  vector  potential  =  Ai^z  such  that  H^  =  V  x  A®,  =  0,  and 

Aj^  =  0  on  the  conductor  and  see  that  A°^  satisfies  Laplace’s  equation  and  the  Dirichlet 
boundary  condition  on  the  conductor  just  as  the  2r— component  of  an  electrostatic  field. 
Therefore  the  potential  for  the  diffracted  magnetostatic  field  satisfies  (4.6)  with  Go 
replaced  by  G%  which  is  the  static  Dirichlet  Green’s  function  for  the  region  y  >  0,  that  is. 


y  >  0-  (4-24) 


The  diffracted  magnetostatic  field,  at  infinity,  behaves  as  if  it  were  scattered  from  a  cylinder 
with  a  finite  cross  section.  This  cissumption  is  confirmed  by  the  fact  that  the  corresponding 
magnetostatic  field  for  the  bump  (see  Chapter  3)  indeed  behaved  as  if  it  was  scattered  by 
a  cylinder  with  finite  cross  section.  From  [20,  sec.  4.4]  it  then  follows  that  in  the  two- 
dimensional  case  A?f  ~  |Vlnr|  =  i  and  ^A?^  ~=  ^  as  r  -+  oo.  Furthermore,  since 
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~  Inr  and  ^G%{r,f’)  ~  ^  as  r  — ►  oo  we  find  that  the  last  integral  in  (4.24)  tends 
to  zero  as  i?  — »  oo.  Hence, 

<(f)  =  jf j,>0  (4.25) 

which  is  the  static  version  of  (4.7).  We  now  use  exactly  the  same  procedure  which  in  the 
time-harmonic  case  led  from  (4.7)  to  (4.16)  and  find  that  (4.25)  leads  to  the  relation 

which  corresponds  to  (4.16).  Similarly  one  finds 

/  Cf},(f,f')KUf')ds'  +  2  /  G»(r, ?')//?,(?-)*'  =  0,  f  6  B  (4.27) 

Jd  Ja 

which  corresponds  to  (4.19).  Because  the  component  of  the  vector  potential  satisfies 

Laplace’s  equation  the  field  /7°  and  the  current  also  satisfy  the  relation  (4.21).  We  have 
now  shown  that  the  magnetostatic  current  and  field  satisfy  the  integral  equations  (4.20), 
(4.21),  and  (4.22)  and  consequently  the  first  term  in  the  low-frequency  expansion  of  the 
current  in  the  dent  equals  this  magnetostatic  current. 

If  we  let  and  be  the  magnetostatic  current  and  field  in  the  magnetostatic  field 
problem  where  the  scatterer  is  situated  in  the  impressed  magnetostatic  field  we  obtain 
from  (4.23)  the  following  low-frequency  expansion  of  the  current  in  the  dent 

K,  =  -2^sin  -f-  Kl  (4.28) 

where  /("J  — >  0  as  Aid  — >  0.  From  the  integral  equations  (4.20),  (4.21),  and  (4.22)  it  follows 
that  7?°^  and  satisfy  the  static  integral  equations 

feA  (4.29) 


I  A'“°(r')dy  +  /  Hf{f)ds'  =  0 

JD  JA 


(4.30) 


£  GUr,r')Kr{r')ds^  ^-2  =  0,  f  e  D.  (4.31) 


(In  [19]  (4.29),  (4.30),  and  (4.31)  are  derived  directly  from  the  Stratton-Chu  formulas.)  In 
summary,  we  have  determined  the  first  term  in  the  low-frequency  expansion  of  the  current 
in  the  dent  and  in  the  next  section  we  integrate  this  low-frequency  current  to  get  the  low- 
frequency  diffracted  far  field. 
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4.1.3  Low-frequency  diffracted  far  field  (TM) 

To  determine  an  expression  for  the  low-frequency  diffracted  far  field  we  will  first  derive  a 
relation  between  the  current  in  the  dent  and  the  diffracted  field  in  F,  the  region  in  the  upper 
half  plane  (y  >  0)  outside  the  image  of  the  dent. 

When  f  e  F  and  f'  G  int{A  U  D)  the  Dirichlet  Green’s  function  satisfies  the 

homogeneous  Helmholtz  equation,  and  Green’s  second  identity  gives 

0  =  f •  (4.32) 

With  the  boundary  conditions  =  0  on  F  and  Gd  =  0  on  4  (4.32)  becomes 

Id  E,{f')—GD{f,f')ds\  f  6  F.  (4.33) 

Use  of  the  expression  (4.7)  for  the  diffracted  field,  the  fact  that  F’  +  F^  =  0  on  4,  and  the 
relation  (4.14)  between  the  current  and  the  electric  field  transform  (4.33)  to 


f  e  F 


(4.34) 


which  is  the  relation  between  the  current  in  the  dent  and  the  diffracted  field  in  F.  This 
expression  was  first  derived  by  As  vestas  and  Kleinman  [18]. 

From  the  asymptotic  expansion  of  the  Hankel  function  one  finds  that 


GD{f,f') 


fcsin^(y' +  0(W)],  r co. 


(4.35) 


The  expansion  (4.35)  and  the  current  expansion  (4.28)  inserted  into  the  field  expression 
(4.34)  gives  the  diffracted  far  field  to  order  (kd)^  as 


Ff(r)  ~  {kdfe  ^  ^2 -^o  sm  ^  sm  </>' 

(4.36) 

where 

Do  =  J  I<f{f)y'ds' 

(4.37) 

which  is  the  final  expression  for  the  TM  low-frequency  diffracted  far  field  of  the  dent  D  in 
the  ground  plane.  We  see  that  this  TM  low-frequency  far  field  is  that  of  a  magnetic  dipole 
in  the  x— direction  [11,  sec.3.8]. 

4.2 

Transverse  electric  (TE)  polarization 

The  incident  electric  and  magnetic  fields  are  here  given  by 

H'{f)  = 

(4.38) 
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and 

E'{f)  =  y^(sin  <l>'x  -  cos  (4,39) 

The  reflected  magnetic  field  is  given  by  (4.38)  with  y  replaced  by  —y.  Similarly,  the  reflected 
electric  field  is  given  by  (4.39)  with  x  and  y  replaced  by  —  i  and  — y,  respectively. 

4.2.1  Integral  equation  for  the  current  (TE) 

We  will  start  by  giving  a  short  derivation  of  Asvestcis  and  Kleinman’s  [18]  coupled  TE  integral 
equations  for  the  current  in  the  dent  and  the  magnetic  field  in  the  aperture  since  we  need 
the  details  of  the  derivation  to  derive  our  low-frequency  results. 

Performing  a  derivation  similar  to  the  one  that  led  to  (4.8)  yields  the  following  relation 
between  the  magnetic  field  in  the  region  y  >  0  and  the  magnetic  field  in  the  aperture. 

H.if)  =  Hi(f)  +  n:(f)  -  so  0.  (4.40) 

Since  the  integral  in  this  equation  is  an  even  function  of  y  we  have 

H,{ri)  =  -h  Hlifi)  -2  f  €  ini{A\jD)  (4.41) 

where  fj  is  the  image  in  the  ground  plane  of  the  point  f  given  in  (4.4).  To  get  rid  of  the 
integral  involving  the  derivative  of  the  magnetic  field  in  the  aperture  we  will  derive  one  more 
relation  containing  this  integral. 

Using  Green’s  second  identity  and  the  fact  that  Hz  satisfies  the  homogeneous  Helmholtz 
equation  we  obtain 

—  f  HJf')-^G{f,f')ds\  r€int(AUD).  (4.42) 

Jaud  on 

Substitution  of  the  integral  in  (4.41)  into  (4.42)  yields 

+Hl{fi)  +  Hl{fi)],  f  G  int{A  U  D).  (4.43) 

Letting  f  approach  A,  using  the  result  of  the  integration  near  the  singularity  from  [14,  App.j 
with  the  boundary  relations 

— G(f ,  f')  =  0,  y  =  y'  =  0,  X  /  x\  (4.44) 
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and 


H,if)  =  Kt{f),  reD 


(4.45) 


(4.43)  converts  to 

This  is  a  relation  between  the  current  in  the  dent  and  the  magnetic  field  in  the  aperture. 
To  determine  the  current  in  the  dent  we  need  one  more  relation  of  this  kind.  Using  Green’s 
second  identity,  the  fact  that  if*  satisfies  the  scalar  Helnaholtz  equation,  and  the  boundary 
condition  that  =  0  on  D  yield  that 

=  H^(f)^GD{f,f*)ds\  f&int{A\jD).  (4.47) 

JAuD  an' 

Letting  f  approach  D  in  (4.47),  using  the  result  of  the  integration  near  the  singularity  from 
[14,  App.],  and  using  (4.45)  gives  us 

iAr.(r)  =  -  £  K,(f')^fio{r,f')ds'-2l^U,{f')^a(f,f')ds',  f  6  D  (4.48) 

which  is  the  second  relation  between  the  current  in  the  dent  and  the  magnetic  field  in  the 
aperture.  The  relations  (4.46)  and  (4.48)  constitute  a  pair  of  coupled  integral  equations  for 
the  current  in  the  dent.  (They  also  can  be  derived  directly  from  the  Stratton-Chu  formulas 
[19].) 

4.2.2  Low-frequency  current  (TE) 

We  will  now  use  the  coupled  integral  equations  (4.46)  and  (4.48)  to  determine  the  first  two 
terms  in  the  low-frequency  expansion  of  the  current  in  the  dent.  Using  the  low-frequency 
expansion  (2.18)  for  -^G  and  the  similar  expansion  for  in  the  integral  equations  (4.46) 
and  (4.48)  we  find  that 

K°{f')-^G°{f,  f')ds’  =  -f  1,  f  €  A  (4.49) 

and 

=  -£  K;{f)-^Gl{f,f')ds'-2  f  e  D  (4.50) 

where  and  are  the  first  terms  in  the  low-frequency  expansions  of  the  current  Kt  and 
magnetic  field  respectively.  In  [15]  it  is  shown  that  these  two  coupled  integral  equations 
have  a  unique  solution. 

We  will  now  show  that  these  are  the  integral  equations  for  the  magnetostatic  current  and 
field  in  the  case  where  the  scatterer  is  situated  in  the  impressed  magnetostatic  field 

HLuc  =  ^’U=o  =  (4.51) 
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To  do  this  we  define  a  reflected  and  a  diffracted  magnetostatic  field  as  in  the  TM  case  (see 
the  discussion  after  (4.23)).  We  find  that  the  reflected  magnetostatic  field  is  z  for  y  >  0 
and  — z  for  y  <  0.  Similarly,  it  is  found  that  the  diffracted  magnetostatic  field  is  2z  in  the 
region  bounded  by  /I  U  /)  and  zero  everywhere  else.  These  observations  make  it  very  easy  to 
use  the  same  procedure  that  gave  us  the  time-harmonic  equations  (4.46)  and  (4.48)  to  show 
that  (4.49)  and  (4.50)  indeed  are  the  equations  for  this  magnetostatic  field  which  simply 
equals  2z  in  the  free-space  region.  Thus  the  current  in  the  dent  and  the  magnetic  field  in 
the  aperture  may  be  expanded  as 


Kt  =  2  +  Kl,  H,  =  2  +  Hl 


(4.52) 


where  and  tends  to  zero  as  fed  —v  0.  We  have  hereby  determined  the  first  term  in  the 
low-frequency  expansion  of  the  current  in  the  dent  and  we  will  now  determine  the  second 
term  in  this  low-frequency  expansion. 

If  we  insert  the  expansions  (4.52)  of  the  current  and  field  as  well  as  the  low-frequency 
expansions  of  the  Green’s  functions  into  the  integral  equations  (4.46)  and  (4.48)  we  get  from 
the  static  integral  equations  (4.49)  and  (4.50) 


f  /\'/(r')^G®(r,r')ds'  =  +jkx  cos  <f>\  f  e  A 

JD  on 


(4.53) 


and 


-  2  f’)ds\  f  6  D  (4.54) 


where  and  Hj  are  the  second  terms  in  the  low-frequency  expansion  of  the  current  in  the 
dent  and  the  magnetic  field  in  the  aperture.  We  therefore  have  the  following  expansions 


Kt  =  2+  kdK°^  -I-  H,  =  2  +  kdH°^  4-  if" 


(4.55) 


where  Kf^  =  and  //“*  =  are  independent  of  kd.  Furthermore,  and 

tend  to  zero  as  kd  — >  0.  Thus,  the  first  two  terms  in  the  low-frequency  expansion  ».f  the 
current  in  the  dent  are  of  order  (A:d)”  and  kd.  Because  the  second  term  in  the  expansion  of 
the  derivative  of  the  Green’s  function  (see  (2.18))  is  of  order  (A:d)^ln  kd  the  third  term  in  the 
expansion  of  the  current  will  in  general  be  a  function  of  both  kd  and  In  kd.  Consequently, 
the  current  does  not  have  a  power  series  expansion. 

Next  we  will  show  that  the  current  A'®*  can  be  found  from  the  solution  to  the  electrostatic 
problem  with  the  scatterer  situated  in  the  impressed  electrostatic  field 

Ktatic  =  cos  <j>'y.  (4.56) 

To  do  this  we  introduce  a  vector  potential  A®  =  F^z  for  the  electrostatic  solution  A®  so 
that 

A®  =  X  (4.57) 
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We  find  that  satisfies  Laplace’s  equation  and  the  Neumann  boundary  condition  on  the 
conductor.  As  we  did  for  the  magnetostatic  problem,  divide  the  scattered  electrostatic  field 
into  reflected  and  diffracted  electrostatic  fields.  The  diffracted  electrostatic  field  at  infinity 
behaves  as  if  it  were  scattered  from  a  cylinder  with  a  finite  cross  section.  This  assumption 
is  confirmed  by  the  corresponding  electrostatic  field  in  the  bump  problem  (see  Chapter  3). 
Because  the  total  charge  on  the  conductor  is  zero  it  follows  from  [20,  sec.3.11]  that  the 
potential  for  the  diffracted  field  satisfies  ~  (Vlnr|  =  |;  and  ~  ^  as  r  — >  oo. 

Using  these  facts  and  the  same  derivation  that  led  to  the  time-harmonic  equations  (4.46)  and 
(4.48)  implies  that  satisfies  these  two  equations  with  Go,  G,  Kt,  Hz,  and  i/*  replaced 
by  G%,  G°,  F^,  F^,  and  F^,  respectively.  Here  Ff'  is  the  potential  for  the  impressed 
electrostatic  field  (4.56).  We  find  that  this  potential  must  be  given  by 

F^'{f)  =  cos  4,'x  +  C  (4.58) 

where  G  is  a  constant.  The  coupled  integral  equations  for  F®  are  thus 

/  F,V)/7<^(FF)ds'  =  -if(f)-|-xcos,^‘-hG,  feA  (4.59) 

Jd  on' 


(4.59) 


(4.60) 


Using  the  divergence  theorem  and  the  result  of  the  integration  near  the  singularity  [14,  App.] 
one  finds 


5^  =  -/. - 2/^cAG»(r-,o*'.  f  €  P. 


(4.62) 


Consequently,  the  constant  C  on  the  right  hand  side  of  the  integral  equation  (4.59)  only 
adds  the  constant  C  to  the  solution  to  (4.59)-(4.60)  with  G  =  0.  Because  a  constant  may 
be  added  to  the  potential  F®  without  changing  the  electrostatic  solution  F®,  let  C  =  0  in 
(4.59).  By  comparing  the  integral  equations  (4.53)-(4.54)  to  (4.59)-(4.60)  with  G  =  0  we 
find  that  Kt^{r)  =  ^F^{r)  when  f  £  D.  This  completes  the  proof  that  the  second  term  in 
the  low-frequency  expansion  of  the  current  can  be  found  from  the  electrostatic  solution. 

It  is  convenient  to  introduce  a  potential  F^^  which  is  the  solution  to  the  equations 

re /I.  (4.63) 

and 

iff-’Cr)  =  -  £  /^0{f')  Ag«  (r,f')*'  -  2  /?°(r')^G»(r,f')Js',  r  e  D  (4.64) 

The  low-frequency  expansion  of  the  current  can  now  be  written  in  terms  of  Ff^  as 

I<t  =  2  +  2jk  cos  <l>'Ff^  -}-  FT  (‘1  -65) 

where  — »  0  as  kd  — >  0. 


(4.64) 


(4.65) 
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4.2.3  Low-frequency  diffracted  far  field  (TE) 


Having  derived  a  low-frequency  expansion  of  the  current  in  the  dent  in  the  previous  section, 
we  will  now  integrate  this  low-frequency  expansion  to  get  the  low-frequency  diffracted  far 
field.  We  start  by  deriving  a  formula  that  presses  the  TE  diffracted  field  in  the  region  F 
in  terms  of  the  current  in  the  dent. 

To  do  this,  note  that  the  Neumann  Green’s  function  G^lf,  f')  satisfies  the  homogenous 
Helmholtz  equation  when  f  E  F  and  f'  E  int{AUD).  Green’s  second  identity  therefore  gives 


0=1  [GN{fy)^H,{f')-H,{f')^GN{f,f')]ds\  fEF.  (4.66) 

JAijD  On  On 

Because  Gn  and  satisfy  the  Neumann  boundary  condition  on  A  and  D,  respectively,  and 
because  Gv  =  2G  on  A,  (4.66)  implies 


j^/^.(r')^G/v(r,r')ds'  =  2jf  G(f,f')^/f.(r0^5',  f  E  F. 
Using  this  in  the  field  expression  (4.40)  gives  [18] 

which  is  the  relation  between  the  diffracted  field  and  the  current  in  the  dent. 
From  the  asymptotic  expansion  of  the  Hankel  function  it  is  found  that 


5  ^  ^ 

— — Gs{fi  r')  ~  A;— 7=— 7=-n'  •  [x  cos  (f)  +  xjkx'  cos^  (f> 
On'  V27r  y/kr 

Fyjky'  s\v^  <j>  +  0{{kdY)\,  r  —*  oo. 


(4.67) 


(4.68) 


(4.69) 


If  the  expansion  (4.65)  of  the  current  and  the  expansion  (4.69)  of  the  Green’s  function  are 
inserted  into  the  field  expression  (4.68)  one  finds  that  to  order  (kd)^  the  diffracted  far  field 
is  given  by 


H^{r)  ~  -(A:d)^^^^^^[^(xx'cos^(^-|-  yy'sin^^6)  •  h'ds 


+  COS  <f> COS  <f>'  f  ^-F^^{r')ds'].  (4.70) 

Jd  Os 


With  the  help  of  the  divergence  theorem  one  finds  that  the  first  integral  in  this  equation 
simply  equals  the  area  of  the  region  bounded  by  AU  D,  that  is,  the  area  of  the  cross  section 
of  the  dent.  The  second  integral  in  (4.70)  may  be  written  as 


Id  =  /„  a'"'  ■  (4.71) 

after  integrating  by  parts  once  and  introducing  the  electrostatic  solution  to  the  problem 
in  which  the  scatterer  is  situated  in  the  impressed  electrostatic  field 
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Before  writing  the  final  expression  for  the  diffracted  far  field,  we  shall  prove  the  relation¬ 
ship 

n-E°^if)  =  Kfif),  fGD  (4.72) 

between  the  TE  electrostatic  solution  in  (4.71)  and  the  TM  magnetostatic  solution  in  (4.28). 

To  do  this  introduce  a  scalar  potential  so  that  =  —Vtp  with  tj;  satisfying  Laplace’s 
equation  and  the  Dirichlet  boundary  condition  on  the  conductor  exactly  as  does  the  vector 
potential  >1°  in  Section  4.1.2.  Furthermore,  because  the  total  charge  on  the  conductor  is 
zero  the  potential  for  the  diffracted  field  also  behaves  like  the  potential  .4^  at  infinity 
[20,  sec.3.11].  Therefore,  the  relations  ^ix  =  and  ^’  =  -y  +  C  where 

C  is  a  constant  and  i/>'  is  the  potential  for  the  impressed  electrostatic  field,  imply  from  the 
coupled  integral  equations  (4.26)-(4.27)  that  ^  satisfies  the  coupled  integral  equations 


and 


r  €  /I, 


=  0,  re  D. 


(4.73) 

(4.74) 


Using  the  fact  that  0  satisfies  Laplace’s  equation  and  employing  the  divergence  theorem  in 
the  region  bounded  by  yl  U  D  give 


/  ■^■ip{f')ds'  =  f  h'  •  V^(f')ds'  =  /  =  0.  (4.75) 

JauD  on'  JaoD  Jint(AuD) 


Using  the  relation  n  •  E^^{f)  =  f  €  D,  and  comparing  (4.73),  (4.74),  and  (4.75) 

to  (4.29),  (4.31),  and  (4.30),  respectively,  shows  that  (4.72)  indeed  is  a  correct  identity. 

According  to  the  relations  (4.71)  and  (4.72),  the  diffracted  far  field  in  (4.70)  may  be 
written  as 

~  +  £>o  cos  <;f> cos  (^']  (4.76) 

where  Ad  is  the  area  of  the  cross  section  of  the  dent  and 


Do  =  (4.77) 

The  magnetostatic  current  is  found  from  the  coupled  integral  equations  (4.29),  (4.30), 
and  (4.31)  and  the  electrostatic  vector  potential  is  found  from  the  coupled  integral 
equations  (4.63)  and  (4.64). 

The  diffracted  far  field  in  (4.76)  is  the  contribution  from  a  magnetic  dipole  in  the 
2— direction  (the  term  with  Ad)  and  an  electric  dipole  in  the  y— direction  (the  cosine  term) 
[11,  sec. 3.8].  Furthermore,  the  constant  in  the  TE  electric  dipole  term  in  (4.76)  is  the  same 
as  in  the  TM  magnetic  dipole  term  (4. .36). 

We  have  now  shown  that  the  low-frequency  scattering  from  a  two-dimensional  dent  in 
a  ground  plane,  as  for  bump  scattering,  reduces  to  the  evaluation  of  a  single  constant  Dq. 
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This  constant  is  determined  from  a  static  solution  and  depends  only  on  the  shape  of  the 
cross  section  of  the  dent. 

Throughout  the  derivation  that  led  to  (4.76)  it  was  assumed  that  the  dent  was  contin¬ 
uously  lined  by  a  conductor  so  that  our  results  in  this  paper  do  not  apply  to  the  slit  in  a 
ground  plane.  However,  since  the  slit  is  the  complement  of  the  strip  it  is  seen  from  Chapter 
2  that  the  low-frequency  fields  diffracted  by  a  slit  for  TE  polarization  must  be  of  order 
where  d  is  the  width  of  the  slit.  Thus,  the  low-frequency  fields  diffracted  by  slits  and  dents 
are  of  different  order  and,  in  particular,  the  slit  scatters  more  strongly  than  the  dent  at 
low  frequencies.  The  physical  reason  for  this  enhanced  scattering  is  that  the  slit  completely 
“stops”  the  currents  on  the  ground  plane  while  the  dent  only  “diverts”  them. 
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Chapter  5 

Evaluation  of  constants  and 
verification  of  low-frequency 
expressions 


This  section  evaluates  numerically  the  constants  occurring  in  the  low-frequency  far-held  ex¬ 
pressions  for  particular  geometries.  Furthermore,  the  expressions  will  be  verihed  by  compar¬ 
ing  them  to  low-frequency  results  obtained  from  exact  time-harmonic  eigenfunction  solutions 
when  such  solutions  exist. 


5.1  Circular  cylinder 


Begin  by  considering  the  circular  cylinder  with  radius  a  situated  with  its  center  at  the  origin 
of  the  coordinate  system.  From  the  exact  eigenfunction  solution  it  is  found  that  the  low- 
frequency  expression  for  the  scattered  far  helds  are  [7,  sec.2.2.1.2,  sec.2.2.2.2] 


TM  : 

Elif) 

TE  : 

Htif) 

lnA:a  y/]^ 

^  2y/^  y/kr 


[tt  -|-  27r  sin  sin  (f>' 

-|-27rcos<^cos<^']. 


(5.1) 


(5.2) 


We  see  that  (5.1)  and  (5.2)  agree  with  the  general  results  in  (2.12)  and  (2.38)  provided  the 
three  constants  in  (2.39)  are 

^  ^  =  27r,  C3  =  0.  (5.3) 

We  can  check  the  expressions  (*’.39)  for  the  three  constants  in  two  ways. 

The  hrst  way  is  to  find  exact  expressions  for  the  electrostatic  charges  and  0“*'  using 
the  method  of  separation  of  variables  and  then  insert  these  into  (2.39).  It  is  ecisily  found 
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that  [21,  pp.  1184- 1185] 


<7°'^  =  2e  cos  =  2e  sin  (j) 


(5.4) 


and  by  inserting  these  into  (2.39)  one  recovers  (5.3). 

The  second  way  of  verifying  the  constants  is  to  solve  the  integral  equations  (2.42)  and 
(2.43)  numerically  and  then  use  (2.47)  to  evaluate  the  constants.  These  two  integral  equa¬ 
tions  were  solved  numerically  by  using  pulse  expajision  functions  and  point  matching.  The 
moment  matrix  elements  were  calculated  using  Simpson’s  rule  of  integration.  The  following 
values  of  the  constants  were  obtained  with  400  pulse  expansion  functions 

^  ^  =  6.283069,  C3  =  0.4  •  10"*^  (5.5) 

which  verify  our  general  low-frequency  results  for  scattering  from  cylinders  with  finite  cross 
section.  We  also  solved  the  integral  equations  (2.40)  and  (2.41)  for  and  and  inserted 
these  numerical  solutions  into  (2.39).  Again  the  result  agreed  with  (5.3). 


5.2  Semi-circular  bump 


Next  consider  the  semi-circular  bump  of  radius  a.  Rayleigh  [17]  used  the  exact  eigenfunction 
solution  for  the  circular  cylinder  to  show  that 


TM  :  £f(f)  ~  (/ca)^e--'^ 

TE  :  H^{f)  ~ 


[2  e"''*’' 

\ - sin  6  sin  6' 

V  TT  y/kr 

/2  TT 

- +  ’rcos<^cos<^’]. 

TT  y/kr  I 


(5.6) 

(5.7) 


We  see  that  this  agrees  with  the  general  low-frequency  expressions  (3.13)  and  (3.19)  provided 
that 


(5.8) 


Again  from  the  method  of  separation  of  variables,  <7®^  =  26  sin  ^  and  one  finds  that  A'®®  = 
— 2sin(^.  Inserting  these  functions  into  (3.14)  one  recovers  the  value  of  Bq  given  in  (5.8). 
One  may  also  combine  the  numerical  results  for  the  circular  cylinder  obtained  by  numerically 
solving  the  integral  equation  (2.42),  with  the  relation  Bq  =  —  |C2  where  C2  is  given  in  (5.3). 
Alternatively,  we  can  solve  the  simple  integral  equation  (3.9)  for  the  current  K°^.  This  was 
done  using  pulse  expansion  functions  and  point  matching.  The  following  value  of  Bq  was 
obtained  with  400  expansion  functions 


^  =  -3.1415965  (5.9) 

which  confirms  our  general  low-frequency  expression  for  the  diffracted  far  fields  of  a  cylin¬ 
drical  bump. 
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5.3  Knife-edge  bump 

Consider  the  vertical  knife-edge  bump  with  height  h.  Using  the  exact  time-harmonic  eigen¬ 
function  solution  for  the  strip  [7,  ch.4]  and  Rayleigh’s  method  of  constructing  an  equivalent 
scattering  problem,  it  is  found  that  the  low-frequency  expressions  for  the  knife-edge  bump 
are 


TM  ; 

TE  : 


Et{r)  ~  (khfe 


H^ir)  ~ 


V  \/^  2 
V  ’T  \/Xr  2 


sin  sin  <t>' 
cos  cos  (f>*. 


(5.10) 

(6.11) 


We  see  that  these  expressions  agree  with  the  general  expressions  (3.13)  and  (3.19)  provided 


that 


-  _5 
“  “2‘ 


(5.12) 


Again  we  shall  verify  that  (3.14)  gives  this  value  for  Bq  in  two  different  ways.  First  we  will 
find  the  static  current  using  separation  of  variables  in  elliptic  coordinates  and  second 
we  will  find  by  solving  the  simple  integral  equation  (3.9)  numerically. 

To  determine  from  the  method  of  separation  of  variables  it  is  convenient  to  use 
the  results  in  Morse  and  Feshbach  [21,  pp. 1195-1197].  is  the  static  current  in  that 
magnetostatic  problem  in  which  a  strip  of  width  2h,  which  has  its  cross  section  directed 
along  the  y-axis,  is  situated  in  the  magnetostatic  field  x.  Therefore,  we  define  the  elliptic 
coordinates  by 

X  = —hsiniifisinO,  y  =  h  cosh  fi  cos  0  (5.13) 


where  /x  >  0  and  0  <  0  <  2ir.  The  ^-coordinate  and  the  strip  given  by  /x  =  0  are  shown  in 
Figure  5.1.  The  scattering  problem  for  can  be  solved  using  a  vector  potential  A°  =  zA'l 
for  the  magnetostatic  field.  The  impressed  magnetostatic  field  can  be  represented  by  the 
impressed  vector  potential 

AP'  —  yz  =  hcosh  HCOS0Z.  (514) 

The  total  vector  potential  A”  must  be  constant  on  the  strip  /x  =  0  and  the  vector  potential 
for  the  scattered  field  must  tend  to  zero  at  infinity.  Therefore  is  given  in  terms  of 
the  elementary  solution  e~^cos0  [21,  p.ll95],  that  is. 


A”*  =  Ce  **  COS0Z 


(5.15) 


where  C  is  a  constant.  Requiring  that  the  vector  potential  for  the  total  magnetostatic 
field  be  zero  on  the  strip  implies  that  C  =  —h,  that  is. 


A°  =  h  sinh  fx  cos  0z 


—X 


COS0  . 

— tZ. 
sm  0 


(5.16) 
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e  =  v 


Figure  5.1:  Strip  of  width  2h  in  an  elliptic  coordinate  system. 


Consequently,  the  current  is 


COS0 

sinS 


0  <  &  <  IT 
TT  <  0  <  27r 


(5.17) 


The  constant  Bq  in  (3.14)  is 

a  =  i  r  yK°^{f)ds  =  -h-‘  r  cos^  M<l  =  (5.18) 

Z  */^=o  ^  Jq  ^ 

which  agrees  with  (5.12). 

Another  verification  is  obtained  by  solving  the  integral  equation  (3.9)  numerically  and 
then  integrating  the  numerical  solution  to  get  Bq.  The  following  value  of  Bo  was  obtained 
with  400  pulse  expansion  functions 

^  =  -1.5706168  (5.19) 

which  again  verifies  the  general  low-frequency  expressions  for  the  diffraction  from  cylindrical 
bumps. 


5.4  Right-angled  cylinders  and  bumps 

We  now  present  numerically  evaluated  constants  for  the  low-frequency  expressions  in  cases 
where  no  exact  time-harmonic  eigenfunction  solutions  exist. 


37 


Using  the  integral  equations  (2.42)  and  (2.43)  with  400  pulse  expansion  functions  it  was 
found  that 

^  =  ^=8.7596,  C3  =  0  (5.20) 

for  a  square  cylinder  with  side  length  2d  and  sides  parallel  to  the  x  and  y  axes.  It  was 
urthermore  found  that  the  constants  changed  to 


^  =  13.3361,  ^  =  25.1541,  ^3  =  0 


(5.21) 


for  a  rectangular  cylinder  with  a  4d  side  length  parallel  to  y  and  a  2d  side  length  parallel  to 

X. 

The  following  bump  constants  were  obtained  by  solving  numerically  the  simple  integral 
equation  (3.9)  using  400  pulse  expansion  functions.  For  the  square  bump  with  side  length 
2d  it  was  found  that 

^  =  -12.5666  (5.22) 


and  for  the  right-angled  triangular  bump  with  height  ^/2d  we  computed 


=  -2.1881684. 


(5.23) 


5.5  Semi-circular  dent 

Consider  the  semi-circular  dent  with  radius  a  situated  with  its  center  at  the  origin  of  the 
coordinate  system.  The  dual-series  eigenfunction  solution  [22]  shows  that  to  order  {kaY 


E^{f)  ~  -{kaYe~^^^*\  —^-y={0.5S)  sin  (j>s\n(l>' 

V  TT  vArr 

Hf{f)  ~  {kaYe~^'^^*\  +  (0.58)  cos  cos 

V  ’T  Vkr  2- 


H\{f)  ~  {kaYe 


These  expressions  agree  with  our  general  result  in  (4.36)  and  (4.76)  provided  that 

^  =  0.58. 


(5.24) 


(5.25) 


(5.26) 


To  determine  the  value  of  Dq  we  solved  numerically  the  coupled  integral  equations  (4.29), 
(4.30),  and  (4.31)  determining  the  magnetostatic  current  K°^.  In  the  numerical  solution  we 
used  pulse  expansion  functions  and  point  matching  and  the  equation  (4.30)  was  incorporated 
by  simply  adding  it  to  the  linear  equations  obtained  from  (4.29)  and  (4.31).  Simpson’s 
rule  of  integration  was  used  to  evaluate  the  moment  matrix  elements.  Using  150  pulse 
expansion  functions  for  both  the  current  in  the  dent  and  the  magnetic  field  in  the  aperture 
(300  unknowns  in  total)  we  found  using  (4.37)  to  determine  Do  that 


^  =  0.58236. 


(5.27) 
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We  also  solved  numerically  the  coupled  integral  equations  (4.63)  and  (4.64)  for  the  elec¬ 
trostatic  vector  potential  .  These  equations  were  also  solved  using  pulse  expansion 
functions  and  point  matching.  Again  Simpson’s  rule  of  integration  was  used  to  evaluate  the 
moment  method  matrix.  Using  225  pulse  expansion  functions  for  both  the  potential  in  the 
dent  and  the  potential  in  the  aperture  (450  unknowns  in  total)  we  found  using  (4.77)  that 

^  =  0.58238.  (5.28) 

a? 

Comparing  the  values  of  Dq  obtained  solving  static  integral  equations  with  the  one  ob¬ 
tained  from  the  dual-series  eigenfunction  solution  [22]  confirms  our  general  low-frequency 
results. 


5.6  Rectangular  dents 

Finally  some  results  will  be  given  for  the  const«int  Dq  in  the  case  of  rectangular  dents. 

The  following  results  were  found  using  the  integral  equations  (4.63)  and  (4.64)  with  225 
pulse  expansion  functions  for  both  the  potential  in  the  dent  and  aperture  (450  unknowns  in 
total).  For  the  rectangular  dent  with  depth  4d,  2d,  d,  0.5d,  and  O.ld  we  found  that  ^  equals 
0.63320,  0.63304,  0.60591,  0.49679,  and  0.16790,  respectively.  The  length  d  in  is  shown  in 
Figure  4.1. 
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Chapter  6 
Conclusion 


We  have  evaluated  the  low-frequency  electromagnetic  scattering  from  the  perfectly  conduct¬ 
ing  cylinder  with  finite  cross  section  and  the  perfectly  conducting  cylindrical  bump  and  dent 
in  a  ground  plane. 

For  the  cylinder  with  finite  cross  section  the  low-frequency  scattered  far  field  for  TM 
polarization  is  independent  of  the  shape  of  the  cross  section  of  the  cylinder  and  of  order 
where  d  is  a  characteristic  dimension  of  the  cylinder.  This  low-frequency  result  is  not  related 
to  a  corresponding  static  field  problem. 

For  TE  polarization  the  scattered  field  is  of  order  {kdy  and  it  consists  of  a  contribution 
from  a  magnetic  dipole  along  the  axis  of  the  cylinder  and  an  electric  dipole  in  a  direction 
normal  to  the  axis  of  the  cylinder.  The  magnetic  dipole  moment  is  found  directly  from  the 
area  of  the  cross  section  of  the  cylinder.  The  electric  dipole  moment  is  found  by  solving  an 
electrostatic  problem,  that  is,  a  two-dimensional  potential  problem,  for  two  impressed  fields 
and  integrating  these  two  electrostatic  solutions  around  the  cylinder.  These  electrostatic 
solutions  are  determined  from  simple  static  integral  equations  and  depend  only  on  the  shape 
of  the  cylinder. 

For  both  the  cylindrical  bump  on  a  ground  plane  and  dent  in  a  ground  plane  the  low- 
frequency  diffracted  field  for  TM  and  TE  polarization  is  of  order  {kd)^,  where  d  is  a  charac¬ 
teristic  dimension  of  the  bump  or  dent.  The  low-frequency  TM  diffracted  far  field  is  that  of 
a  magnetic  dipole  normal  to  the  axial  direction  and  parallel  to  the  ground  plane. 

The  low-frequency  TE  diffracted  far  field  for  both  the  cylindrical  bump  and  dent  consists 
of  a  contribution  from  a  magnetic  dipole  in  the  axial  direction  and  an  electric  dipole  normal 
to  both  the  axial  direction  and  the  ground  plane.  The  TE  magnetic  dipole  moment  is  found 
directly  from  the  area  of  the  cross  section  of  the  bump  or  dent. 

Both  the  TM  magnetic  dipole  field  and  the  TE  electric  dipole  field  can  be  written  as 
a  constant  times  a  known  simple  function.  It  is  proven  that,  remarkably,  this  constant, 
which  depends  only  on  the  shape  of  the  bump  or  dent,  is  the  same  for  both  the  TM  and  TE 
polarizations.  This  constant  can  be  found  by  solving  either  a  magnetostatic  or  electrostatic 
field  problem,  that  is,  a  two-dimensional  potential  problem,  and  performing  an  integration  of 
these  solutions  over  the  bump  or  dent.  This  means  that  both  the  TM  and  TE  low-frequency 


40 


clifTracted  far  fields  for  an  arbitrarily  shaped  bump  or  dent  are  completely  determined  by 
calculating  a  single  constant  for  the  bump  or  dent. 

The  low-frequency  expressions  were  confirmed  from  exact  time-harmonic  eigenfunction 
solutions,  and  constants  were  evaluated  for  a  number  of  geometries.  As  mentioned  in  the 
Introduction,  the  low-frequency  expressions  derived  in  this  report  are  given  in  closed  form 
and  can  therefore  be  used  directly  to  determine  incremental  length  diffraction  coefficients 
for  calculating  the  scattered  fields  from  curved  narrow  ridges  and  channels  in  conductors 
[9].  These  incremental  length  diffraction  coeflBcients  determined  from  the  low-frequency 
expressions  derived  in  this  report  have  been  applied  recently  to  calculate  the  effects  of  ridges 
and  channels  on  the  fields  of  a  reflector  antenna  [23]. 
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Appendix  A 

Existence  of  an  electrostatic  vector 
potential 


Let  be  the  total  electrostatic  field  in  the  TE  electrostatic  problem  in  which  the  scatterer 
is  the  cylinder  with  finite  cross  section  shown  in  Figure  2.1.  From  Maxwell’s  equations  we 
find  that  has  zero  divergence,  that  is, 


V^  =  0 


(A.1) 


and  because  the  total  charge  on  the  conductor  is  zero  we  have 

j^E°-nds  =  Q  (A.2) 

where  S  is  the  curve  that  describes  the  cross  section  of  the  cylinder. 

We  will  show  that  (A.l)  and  (A.2)  imply  that  there  exists  an  electrostatic  vector  potential 


P  =  F^z 


(A.3) 


so  that 


^  =  V  X  F®. 


(A.4) 


We  start  by  proving  that 

f  EP  •  hds  =  0  (A. 5) 

J  c 

for  every  closed  curve  C  in  the  free-space  region  of  the  x  —  y  plane  of  Figure  2.1.  Here 
h  =  t  X  z  where  t  is  the  tangent  unit  vector  to  C. 

Consider  first  the  case  where  C  does  not  enclose  the  cylinder  S:  Then,  since  the  electro¬ 
static  field  has  zero  divergence  it  follows  from  the  divergence  theorem  that  (A.5)  is  satisfied. 

Consider  then  the  case  where  C  encloses  5;  The  divergence  theorem  shows  that  the 
integral  over  C  equals  plus  or  minus  the  integral  over  5,  which  by  (A.2)  is  zero. 

We  have  now  shown  that  the  integral  over  every  closed  free-space  curve  C  of  n  •  F®  is 
zero. 
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Let  {E°,Ey,0)  be  the  rectangular  components  of  the  electrostatic  field  E°,  then  we  have 
the  following  identity  on  C 

h  E^  =  i-  [-E^x  +  E^y].  (A.6) 

Because  of  (A. 5)  we  know  that  the  field  —E^x  +  is  conservative  and  it  follows  then  from 
Phillips  [24,  p.52]  that  there  exists  a  scalar  potential  0  such  that 

-  E°x  +  Ely  =  VV>.  (A.7) 

By  inspection  one  finds  that  (A. 4)  holds  with 

#?  =  V-  (A.8) 

and  the  existence  of  a  vector  potential  for  the  electrostatic  field  satisfying  (A.l)  and 
(A.2),  is  proved. 
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Appendix  B 

Dipole  moment  reciprocity  theorem 


Consider  the  two  electrostatic  problems  where  the  cylinder  in  Figure  2.1  is  situated  in  the 
impressed  electrostatic  fields  E{  =  x  and  E2  =  y.  Let  the  solution  to  these  two  static 
problems  be 

=  +  (B-l) 

where  E{  2  are  the  scattered  fields. 

We  will  prove  the  dipole  moment  reciprocity  theorem: 


j  yh  ■  E^ds  =  J^xh-  E^ds 


(B.2) 


which  states  that  the  y  component  of  the  dipole  moment  in  tL„  scattering  problem  with 
E'  =  X  equals  the  x  component  of  the  dipole  moment  in  the  scattering  problem  with  E'  =  y. 
To  prove  this  we  introduce  the  scalar  potentials  ^1,2  such  that 

{B.3) 

The  potentials  V’l  2  for  the  impressed  fields  are  given  by 

=  -X  +  Cl  (B.4) 

and 

,/;«=_y  +  C2  (B.5) 

where  Ci  and  C2  are  constants.  Since  the  total  charge  on  the  conductor  is  zero,  that  is, 


J  n  E^,,ds  =  0 


(B.6) 


it  is  seen  that  (B.2)  is  equivalent  to 


(B.7) 


47 


Because  the  potentials  0}  2  satisfy  Laplace’s  equation  inside  S  we  obtain  from  Green’s  second 
identity 


(B.8) 


If  we  require  that  0i_2  =  0  on  5  we  see  from  (B.8)  that  (B.7)  and  therefore  also  (B.2)  is 
equivalent  to 


(B.9) 


where  0J  2  =  0i,2  ~  01,2  potentials  for  the  scattered  electrostatic  fields. 

Let  Sft  be  the  circle  with  radius  R  and  center  at  x  =  y  =  0.  Assume  that  encloses 
S.  Then  since  0J  2  satisfy  Laplace’s  equation  outside  S,  we  get  from  Green’s  second  identity 
that 
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ds. 


(B.IO) 


Now,  the  total  charge  is  zero,  so  it  follows  from  [20,  sec.3.11]  that  for  two-dimensional 
problems 
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(B.ll) 

(B.12) 


The  asymptotic  conditions  (B.ll)  and  (B.12)  imply  that  the  last  integral  in  (B.IO)  goes  to 
zero  as  R  goes  to  infinity.  Thus,  (B.9)  and  therefore  the  dipole  moment  reciprocity  theorem 
is  proven. 
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